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The objective of the present research is to develop improved turbulence models for the 
computation of complex flows through turbomachinery passages, including the effects of 
streamline curvature, heat transfer and secondary flows. 

Advanced turbulence models are crucial for accurate prediction of rocket engine flows, 
due to existance of very large extra strain rates, such as strong streamline curvature. Numerical 
simulation of the turbulent flows in strongly curved ducts, including two 180-deg ducts, one 
90-deg duct and a strongly concave curved turbulent boundary layer have been carried out with 
Reynolds stress models (RSM) and algebraic Reynolds stress models (ARSM). The RSM & 
ARSM models are successful in the prediction of damping effects of convex curvature. 
However, both models underpredict the turbulence amplification caused by strong concave 
curvature. In order to capture this amplification of turbulence, the time scale (for spectral 
energy transfer) in the dissipation rate (e) equation must be modified. A detailed analysis has 
been carried out for the modifications to the e-equation. An improved near-wall pressure-strain 
correleration has been developed for capturing the anisotropy of turbulence in the concave 
region. 


A comparative study of two modes of transition in gas turbine, the by-pass transition 
and the separation-induced transition, has been carried out with several representative low- 
Reynolds-number (LRN) k-e models. Effects of blade surface pressure gradient, freestream 
turbulence and Reynolds number on the blade boundary layer development, and particularly the 
inception of transition are examined in detail. The present study indicates that the turbine blade 
transition, in the presence of high freestream turbulence, is predicted well with LRN k-e 
models employed. 

The three-dimensional Navier-Stokes procedure developed by the present authors has 
been used to compute the three-dimensional viscous flow through the turbine nozzle passage of 
a single stage turbine. A low Reynolds number k-e model and a zonal k-e/ARSM (algebraic 
Reynolds stress model) are utilized for turbulence closure. The algebraic Reynolds stress 
model is used only in the endwall region to represent the anisotropy of turbulence. For the 
turbine nozzle flow, comprehensive comparisons between the predictions and the experimental 
data obtained at Penn State show that most features of the vortex-dominated endwall flow, as 
well as nozzle wake structure, have been captured well by the numerical procedure. An 
assessment of the performance of the turbulence models has been carried out. The two models 
are found to provide similar predictions for the mean flow parameters, although slight 
improvement in the prediction of some secondary flow quantities has been obtained by the 
ARSM model. It's found that the wake profiles inside the endwall boundary layers arepredicted 
better than those near the mid-span. 
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Objective: 


To develop turbulence models for prediction of turbine flow & thermal 
fields including effects of curvature, rotation and high temperature 


Outline: 


• Introduction 

• Numerical Technique & Turbulence Models 

• 3-D Navier-Stokes Comp, of turbine nozzle flow 

• Turbulence Modeling for Strongly Curved Shear Flows 

• Computation of Turbine Blade Transition and Heat Transfer 

• Concluding remarks 


3-D NAVIER-STOKES PROCEDURE 


Explicit 4-Stage Runge-Kutta Scheme 

Central differencing + smoothing (eigenvalue & local vel. scaling) 
Turbulence Models: 

- Differential Reynolds stress model (high Re no. & low Re no.) 

- Algebraic Reynolds stress model 

- Nonlinear k-e model 

- Two eq. models (low & high-Re-no. versions) 

- e-modification for strong streamline curvature 

Boundary Conditions 

- Characteristic boundary conditions 

- Quasi-3D non-reflecting boundary conditions 

Acceleration Schemes 

- Local time stepping, implicit residual smoothing 

- Implicit treatment of k-e equations 

- Multigrid 2 



Strongly curved shear flows investigated: 


Flow 

Author 

Re 

8/R 

Concave TBL 

Barlow & Johnston 

3.3x104 

0.06 

90-deo duct 

Kim & Patel 


0.05/0.04 

180-dea duct 

Monson et al. 


0.7/0. 2 

180-dea duct 

Sandborn 


1.0/0. 3 

180-dea duct 

Monson et al. 

1.0x10® 

0.7/0. 2 
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Modeling for Curved Shear Flows 


• Modeling of curved flows (mostly mild curvature, convex curvature) 

U/r 

- Mixing-length model: Prandtl's hypothesis F = 1 -a gu / g p 

£ 2 

- k-e model: Launder et al (1977) -C £ 2 (1- C c Rij) , etc. 

- RSM & ARSM: Irvin & Arnot Smith 1975, Gibson & Rodi 1981, etc. 
(used in boundary layer codes, mild/convex curved flows) 

• Comps, of strongly curved flows, e.g., 180-deg duct flows. 

- Monson et al. (1990); 

- Avva et al. (1990); 

- Shih et al. (1994), etc. 

Agreement not satisfactory 

=> Further modeling work on strongly curved flows 



Modifications to e-equation 


Standard e equation: 
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standard values: C^=0.09, C e i=1.44, C e 2=1.92, ak=1.0, o e =1.3 


• Modification in the sink term (Launder et ai. 1977) 

i 

C e2 =C £ 2(1-0.2Rit) 


where Ri 
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Modification of time-scale in the source term (Lumley 1992) 
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Nonliear k-e model (Shih, Zhu & Lumley 1993) 


u.u.=^k5..-v f (U. .+U..) 
i J 3 ij I i,j J,i 
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ARSM (ALGEBRAIC REYNOLDS STRESS 

MODELS) 


• Reynolds stress transport eq.: 

^ij “ ®ij “ *\j + ^ij “ E ij 

• ARSM assumption (Rodi, 1976): 


• Present ARSM (Derived from Gibson 
& Launder RSM (1978) for compressible 
flows) 


ti ff 


P u i u j=-pk 


(P.. -2P5../3X1- C,) +()).. 

U iJ/ v 2' y ij,w 


P+petCj-Dl-ISypk 


p ij = — p u j' u ^3uj/3x k — P u j u ^3uj/dxk and P = P../2 
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Low Revnolds-number RSM Models 


Shima. 1988 

• Based on LRR RSM 

• conventional damping function f w ( f w = exp(- 
(0.015k 1/2 y/n) 4 ) 

• This LRN model reduces to its high-Re version (i.e., LRR model) 
away from the wall. 

Launder & Shima. 1989 

• Based on LRR-Gibson-Launder RSM 

• Use independent Reynolds stress invariants 

• Constants c,, C 2 , C, and c' 2 are functions of the turbulence 
anisotropy parameters 

•This model may not reduce to its high Re version away from the 
wall 



REYNOLDS-STRESS MODEL (RSM) AND 
ALGEBRAIC REYNOLDS-STRESS MODEL (ARSM) 


• Reynolds stresss transport equation : 


8u.u. 

U — U = - 


=-u i u^ jik -ii ] n^ a+ £(„ IJ+ Uj,i)-5- 
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• Models employed in present computations: 

• RSM model = LRR Model (Launder, Reece & Rodi, 1975) 
with Shima near-wall low-Reynolds-number functions 

• ARSM model = Algebraic form of LRR model with Gibson- 
Launder near-wall Pressure-strain correlation 
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# of iter vs. log(res) 


EULER SOLUTION 


PERDICHIZZI TURBINE 

NOZZLE CASCADE 

Chord length 

55.2 mm 

Axial length 

34.0 mm 

Aspect ratio 

1.47 

Inlet blade angle 

76.1 deg 

Outlet blade angle 

14.5 deg 

Inlet Mach number 

0.15 

Outlet Mach number 

0.70 

Reynolds number 

0.84x106 



Nb 


ORIGINAL 






MULTIGRID 


IMPROVEMENT 


5 : 1 Convergence Rate 
4.3 : 1 CPU Time 



iter# 


24 Apr 95 1 1:34:01 
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Tu 


6 . 0 


2.0 


y 


R.Q 



0.0650 

0.0600 

0.0550 

0.0500 

0.0450 

0.0400 

0.0350 

0.0300 

0.0250 

0.0200 

0.0150 

0.0100 


4.0 


0.0 


4.0 


8.0 


23 Apr 95 20:16:0 





440 



0,0650 

0.0600 

0.0550 

0.0500 

0.0450 

0,0400 

0.0350 

0.0300 

0.0250 

0.0200 

0.0150 

0,0100 


23 Apt 95 19:55:08 










0.5{<uu>+<vv>)/U m xl 000 0.5(uu+vv)/U m * 1 000 

6=30 deg 0=90 deg G=180 deg 

Fig. Variation of 0.5(<uu>+<vv>) profile of 180-deg duct flow 








Fig. Profiles at 0=J 


1 . 0 - 



0 . 4 - 


O Data8/R = 50% 

5/ R = 1 5 % 

=25 

=50 


0=90 deg 


9=180 deg 


Fig. Effects of 5/R on mean velocity profile of 180-deg duct flow 
(Simulation by RSM model, data by Sandborn) 


0.5(<UU>+<vv>)/U m xl 000 

x= -1.73 H 0=90 deg 0=180 deg 

Fig. Effects of 8/R on turbulence energy profile of 180-deg duct flow 
(Simulation by RSM model, data by Sandborn) 



Conclusions on Turbulence Modeling for Strongly Curved Flows 

• RSM model provide best predictions for major features of the highly 
curved duct flows, including major attenuation of turbulence near the 
convex wall, strong enhancement of turbulence near the concave wall 
and the extensive separation downstream of the bend. 

• Modeling convex curvature effects is different from modeling concave 
curvature effects, even qualitatively. RSM model is very successful in 
modeling convex curvature, while still underpredicts concave curvature 
effect. 

• Turbulence damping due to convex curvature are also captured well 
by ARSM & Nonlinear k-e model. The isotropic k-e model fails to account 
for this effect and underpredicts the extent of separation. 

• ARSM model is superior to nonlinear k-e model for the curved duct 
flows investigated. The nonlinear k-e model does not capture any 
turbulence enhancement in concave region. 


• All the models provide too slow a recovery process from separation 
downstream of the bend, indicating defects in the modeling of turbulent 
diffusion as well as dissipation terms. 

• Simulation studies & data indicate that the flow inside the bend is not 
sensitive to the upstream inflow conditions, different 5/R leading to only 
minor variation in downstream velocity profiles. 

• To capture concave curvature effect, e-eq. must be modified: 

- Model of Launder et al. & Lumley provide some improvement. 

- Further improvement of e-eq. is needed. 


Design Features of Penn State Turbine Nozzle 

Hub tip ratio 0.7269 

Tip radius 0.4582 m 

Chord(tip) 0.1768 m 

Spacing(tip) 0.1308 m 

Turning angle 70 deg 

Vane Re(outlet) (9-10)x105 
Exit Mach 0.27 



Fig. 3a Computational grid for PSU turbine nozzle 
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Fig. 3b Measurement locations for the turbine nozzle 

(FHP: Five hole probe; 

LDV: Laser doppler velocimeter) 





Fig. Pitchwise-mass-averaged parameters 
(tangential velocity and secondary kinetic energy) 
predicted using 121x73x89 and 99x69x69 grids with 

k-c model 
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Mg 


Computation by k-e/ARSM model 


Measurement 


Fig. 7 Nozzle hub wall static pressure coeffi. (C p ) 





Fig. Pitchwise-averaged total velocity profile 
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Pitchwise-averaged total pressure loss coeff. 
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Fig. 18b Axial ve 


Conclusions 


• Most features of the vortex-dominated endwall flow in the annular 
turbine nozzle have been captured accurately by the 3-D Navier-Stokes 
prediction. The passage-averaged properties, particularly the yaw angle 
and velocity profiles, are captured very well by the present numerical 
computation. 

• The predictions by the anisotropic ARSM model are close to those 
by the isotropic k-e model for the mean flow properties, although slight 
improvement in the prediction of secondary flow (e.g., the secondary 
kinetic energy) has been obtained by the ARSM model. 

• The turbine nozzle secondary flows are primarily driven by pressure 
gradients. The anisotropy of turbulence becomes important when the 
secondary flow rolls up into a distinct vortex. Its dissipation and diffusion 
may only be captured by the ARSM and other anisotropic turbulence 
models. 


• The wake profiles inside the endwall boundary layers are predicted 
better than those near the mid-span. The width and depth of the wake at 
the mid-span are overpredicted due to a premature transition predicted 

by the k-e model on the blade suction surface in the presence of low 
freestream turbulence. The discrepancy in the wake profile are also due 
to the downstream rotor influence. 


oo 



PSU Rotor (midspan) 

• Axial chord = 9.114 cm 

• True chord = 11.13 cm 

• Flow turning angle = 110 deg 

• Re (exit) = 5-7x10^ 

• Mach (exit) = 0.27 
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co 1.0 


Fig. 

(Mea 



Utot/Um 

/ A' 



(computation with k-£/ARSM model) 


flowfield 



.3H 


. 2 - 

y (pitchwise) 

x/C=102% x/C=112% 

Fig. Wake profile at midspan for PSU rotor flow 
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ries of Mark II and C3X: 


Mark II C3X 

63.69 59.89 

70.96 72.38 

12.97 11.77 

13.62 14.49 
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Analysis of Heat Transfer for Mark II & C3X Turbine Nozzle Guide Vanes 
with conditions (Re, Tu, TO, Tw) close to real engine conds. 

• Preds. with engr. accuracy obtained by CH, LB, FL k-e models 

- LB prediction appears to be the best 

- CH performs well in fully turbulent region, 
but tends to smear out transition process, 

also not good for separation-induced transition 

- FL yield delayed transitions for accelerating turbine flows 

• Separated-flow transition lead to much sharper increase of 
heat transfer than the nominal by-pass transition 

• With minimum smoothing and good LRN k-e model, 

2-D N-S method provide good pred. of blade boundary layer 
development, transition & heat transfer under diff. Re, Tu, etc. 
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Future Efforts 


• Turbulence modeling: 

- Investigation of combined effects of curvature & rotation on 
turbulent flowfield in 3-D rotor flows 

- Modeling the source term in e-equation to capture strong concave 

curvature; couple this with RSM 

- Modeling the e-eq. & turbulent diffusion to improve the prediction 
of recovery process after re-attachment 

• Code Development: 

- Multigrid solution of Reynolds stress transport equations 

- Development & implementation of full 3-D non-reflecting B.C. 

• Validation & Simulation: 

- Modeling strong concave curved TBL & duct flows 

- Navier-Stokes simulation of rotor flow with non-reflecting 
boundary conditions & advanced turbulence models 



